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ABSTRACT 

We present an analysis of the distribution of Ha emission measures for the warm ionized medium 
(WIM) of the Galaxy using data from the Wisconsin H- Alpha Mapper (WHAM) Northern Sky Survey. 
Our sample is restricted to Galactic latitudes |6| > 10°. We removed sightlines intersecting nineteen 
high-latititude classical H II regions, leaving only sightlines that sample the diffuse WIM. The distri- 
bution of EM sin \b\ for the full sample is poorly characterized by a single normal distribution, but is 
extraordinarily well fit by a lognormal distribution, with (log EM sin |6|( pc cm -6 ) -1 ) = 0.146± 0.001 
and standard deviation <7 logElvlB m |6j = 0.190 ± 0.001. (logEMsin|6|) drops from 0.260 ± 0.002 at 
Galactic latitude 10 < \b\ < 30 to 0.038 ± 0.002 at Galactic latitude 60 < \b\ < 90. The distribution 
may widen slightly at low Galactic latitude. We compare the observed EM distribution function to the 
predictions of three-dimensional magnetohydrodynamic simulations of isothermal turbulence within 
a non-stratified interstellar medium. We find that the distribution of EMsin|6| is well described by 
models of mildy supersonic turbulence with a sonic Mach number of ~ 1.4 — 2.4. The distribution 
is weakly sensitive to the magnetic field strength. The model also successfully predicts the distribu- 
tion of dispersion measures of pulsars and Ha line profiles. In the best fitting model, the turbulent 
WIM occupies a vertical path length of 400 — 500 pc within the 1.0 — 1.8 kpc scale height of the 
layer. The WIM gas has a lognormal distribution of densities with a most probable electron density 
n p k ~ 0.03 cm" 3 . We also discuss the implications of these results for interpreting the filling factor, 
the power requirement, and the magnetic field of the WIM. 
Subject headings: ISM: structure — turbulence — MHD 



1. INTRODUCTION 

The warm ionized medium (WIM, sometimes called 
the diffuse ionized gas, or DIG) is a major component of 
our Galaxy, consisting of a pervasive, diffuse plasma layer 
with temperatures near 8000 K, a scale height of about 
1 kpc (more tha n three times t he scale height of the 
neutral hydrogen; lRevnoldsHl989[). and a spac e-averaged 
midplane density of 0.03 cm" 3 (|Ferrie rc 2001). By com- 
binin g pulsar dispersio n measures with Ha emission mea- 
sures. iRevnoldsl (|199lD characterized the gas with a sim- 
ple model consisting of discrete clumps with densities 
of 0.08 cm -3 occupying ~ 20 — 40% of the volume of 
a 2 kpc thick, plane parallel layer about the midplane 
of the Galaxy. The WIM is distinct from a component 
of warm, ionized gas associated primarily with classical 
H II regions, which has a scale height of 40 — 70 pc and a 
space-averaged density comparabl e to that of the WIM 
(Fcrriere 2001: iGomez et al"]l2001l ). Thus, in directions 
away from the Galactic midplane, emission from warm, 
ionized interstellar gas is dominated by the WIM. 

The ionization power requirement of the WIM is 
2n 2 a^fhwiM ~ 4 x 10 6 ionizing photons s" 1 cm" 2 or 
~2x 10~ 4 erg s" 1 cm" 2 in the Galactic disk ([Reynolds! 
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1990), where is the recombination coefficient to 

levels > 2, and n, f and /iwim are the electron den- 
sity, volume filling fraction, and scale height of the 
WIM, respectively. With a core colla pse supernova rate 
of one 10 51 erg event per 50 years (jDiehl et al.l l2006h 
in the Galactic disk (r ~ 15 kpc), supernovae inject 
barely enough energy to ionize the WIM even if all 
of the energy from supernovae goes into ionizing the 
WIM. Lyman continuum radiation from O stars injects 
3 x 10 7 ionizing photons s" 1 cm" 2 (|AbbottHl982T) : this is 
the only known source that injects ample energy to ionize 
the WIM. With the O stars and supernovae concentrated 
in the plane, the mechanism by which the ionizing radi- 
ation and shocks reach the 1 kpc scale height is not well 
understood. 

Given the observed velocities and expected sources of 
viscosity , one expects th e WIM to have a high Reynolds 
number (|Beniaminlll999f ) , meaning the gas should be tur- 
bulent. At sm all scales, observa t ions h ave confirmed this 
expectation. lArmstrong et alj (|1995f ) used radio mea- 
surements of interstellar scintillation and fluctuations in 
dispersion measure and rotation measure to demonstrate 
that the ISM is a turbulent medium with density fluctu- 
ations following a Kolmogorov-like power law over scales 
from 2 x 10 s cm (10~ 5 AU) to at least 10 15 cm (70 AU). 
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Their data are consistent with a power law spectrum to 
scales as large as 10 20 cm (30 pc). 

Astrophysical turbulence does not nee d to be 
Kolmogorov-like, however. (See review by iLazarianl 
2006.) Integrated intensity fl uctuations of 13 CO and 
12 CO show a shallow spectrum (jStutzki et al.lll998l ). Us- 
ing models of radiative tran sfer in a turbulent medium, 
lLazarian fc Pogosvanl (|2004D interpreted these fluctua- 
tions as corresponding to a spectrum of turbulent den- 
sity E{k) r-j k~P with a spectral index (3 = 0.8, whereas 
= 5/3 for Kolmogorov turbulence. I Swift fc Welch! 
(2008) found (3 = 1 — 2 for fluctuations in the optically 
thin isotope C 18 in t he low mass, star-formin g molec- 
ular cloud Lynds 1551. IChepurnov et al.l (|2006l ) showed 
that H I data for high latitudes correspond to a velocity 
spectrum steeper than the Kolmogorov value and a den- 
sity spectrum shallower than Kolmogorov. This is well 
in agreement with magnetohydrodynamic (MHD) nu- 
merical simulations (|Beresnvak et al.l l2005| ) which show 
that supersonic turbulence te nds to have s teep v elocity 
and shallow density spectra. iKim fc Rvul (|2005l) found 
similar shallow density spectra at high Mach number in 
hydrodynamic simulations and estimated a sonic Mach 
number of order unity for the WIM. 

Pulsar dispersion measu res provide the column den - 
sity through the g a s, whi ch iManc hestcr fc Tavlorl (1 19811) , 
iTavlor fc CordesT (fl993l ). and iGolnez et alJ (j2001h . 
among others, used to construct models of the distribu- 
tion of warm, ionized gas in the Galaxy. The bulk of the 
information about the temperature, ionization state, and 
kinematics of the WIM has come from the study of faint, 
optical emission lines such as Ha (6563 A), [S II] A6716, 
and [N II] A6583 with Fabry-Perot spectrometers such 
as the Wisconsin H- Alpha Mapper (WHAM). Velocity- 
resolved observations of Ha emission reveal the distri- 
bution and kinematics of the gas, while the line ratios 
[N II] /Ha and [S II] /[N II] are sensitive pri marily to tem- 
perat ure and ionization state, respectively ijHaffner et al.l 
1999). The comparison of Ha emission measure with 
pulsar dispersion measures along the same si ghtline al- 
low e stimates of the "dumpiness" of the gas (| Reynolds! 

EMU). 

In this paper, we show that the Ha emission measure 
of the WIM is extremely well characterized by a lognor- 
mal distribution, demonstrate that the observed distri- 
bution is plausibly produced by MHD turbulence in an 
approximately isothermal medium, and pursue observa- 
tional and theoretical ramifications of this result. In § [2l 
we discuss the WHAM data used in this work and our 
masking of the WHAM data to study only the diffuse 
WIM. In § [3l we demonstrate that the diffuse WIM has 
a lognormal distribution of emission measure, and we 
discuss the occurrence of lognormal distributions in con- 
nection with turbulence in § |4] We compare the results 
of isothermal MHD models of turbulence with the data 
in § [5j finding a reasonable correpsondance. In § [6l we 
use these models to make and then test additional pre- 
dictions for the WIM observations. Finally, in § [7j we 
explore the ramifications of this work on the filling frac- 
tion, power requirement, and magnetic field structure of 
the WIM. 



2. WHAM DATA 



TABLE 1 

Large Classical H II regions, 5 > -30°, |£>| > 10° 





Ionizing Star 






Radius 




Spectral 






Removed 


Name 


Type 


I (°) 


6(°) 


(°) 






Classical H II 


regions 






COph 


09 V 


6.3 


23.6 


7 


HD 161056 


B1.5 V 


18.7 


11.6 


2 


7r Aqr 


Bl Ve 


66.0 


-44.7 


3 


10 Lac 


09 V 


96.6 


-17.0 


1 


AO Cas 


09 IIInn+... 


117.6 


-11.1 


3 


PHL 6783 a 


sd:B 


123.6 


-73.5 


3 


<fi Per 


B2 Vpe 


131.3 


-11.3 


8 


a Cam 


09.5 Iae 


144.1 


14.0 


6 


£ Per 


07.5 Hie 


160.4 


-13.1 


10 


HD 41161 


08 V 


165.0 


12.9 


5 


A Ori (HD 36861) 


08 III 


195.0 


-12.0 


7 


<t Ori 


09.5 V 


206.8 


-17.3 


9 


PC 1034+001 


DO 


247.6 


47.8 


5 


PG 1047+003 


sd:B 


250.9 


50.2 


5 


Spica (a Vir) 


Bl III-IV+... 


316.1 


50.8 


12 


7T Sco 


Bl V+... 


347.2 


20.2 


2 


S Sco 


B0.2 IVe 


350.1 


22.5 


9 




Other masked 


regions 






Gum Nebula 




255.8 


-9.1 


~ 40 


Orion-Eridanus bubble 


45. 


-25. 


~ 25 



Note. — Ionizing stars of H II regions which we removed 
from the WHAM Northern Sky Survey to leave a sample of 
only the diffuse WIM. This list includes only high-latitude 
(|6| > 10°) regions. Some sightlines are contained in mul- 
tiple, overlapping H II regions. The radius removed is the 
maximum extent of significantly enhanced Ha and may be 
larger than the size of asymmetrical H II regions. Coordi- 
nates and spectral types are from the SIMBAD Astronomical 
Database (http : //s imbad.u-strasbg . f r/ ); see also Table 1 of 
IMadse n et al. (200(j). 

a Sec Haffncr (2001) for a discussion of the ionizing sources in 
this region. 

The Wisc onsin H-Alpha Map per (WHAM) Northern 
Sky Survey (IHaffner et al.l l2003) is a spectroscopical map 
of Ha emission in the northern sky within ±100 km s _1 
of the local standard of rest with an angular resolution 
of 1°. The spectra were obtained using a dual-etalon 
Fabry-Perot spectrometer with a spectral resolution of 
12 km s _1 and a sensitivity of 0.15 R (0.34 pc cm -6 ). 
The spectral resolution of the WHAM survey allows the 
subtraction of numerous atmospheric emission lines, in- 
cluding the relatively bright geocoronal Ha line, to estab- 
lish a stable zero point for the interstellar Ha intensities. 

We calculated emission measures from the WHAM sur- 
vey irrtegn^doyer velocities \v\ < 100 km s -1 according 
to (|Revnoldslll9"9l[) 

where n e is the local electron density, In a is the Ha in- 
tensity (1 R = 10 6 (47r) _1 photons cm~ 2 s _1 sr^ 1 ), and 
T is the temperature. To determine the temperature, 
we consider the ionization states in the WIM. Through 
the WIM, [S II] /[N II] line ratios are relatively constant, 
which implies a uniform ionization state, but [S II] /Ha 
and [N II] /Ha line ratios tend to b e higher in lower Ha 
intensity regions and at higher \z\ (|Haffner et al.|[l999l) . 
In combination, these observations imply a range of tern- 
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Fig. 1. — Vertical component of Ho emission measure EM sin |fe| 
a) and with sightlincs intersecting classical H II regions removed 
increments in latitude and 30° increments in longitude. Single p: 
removed in the processing of the survey. 

og EM sin Ibl (pc cm~ 6 ) 



-0.80 -0.30 0.20 0.70 






Fig. 2.— Map of EMsin|6| from the WHAM survey for the 
northern polar region. Rectangular regions with discontinuities 
in intensity represent 7x7° observational "blocks" which were 
observed together; a t the low intensities observed at high latitudes 
(Haffncr et al. 2003), block-to-block changes in the sky subtraction 
represent a significant fraction of the emission. 

peratures in the WIM from 6000 — 10000 K, with higher 
temperatures higher above the plane and in low density 
regions. We adopt a temperature of 8000 K. 

Within about 5° of the midplane, extinction by dust 
obscures many sightlines. Therefore, we restrict our 
data set to sightlines with Galactic latitude |&| > 10° , 
where extinction is negligi ble at Ha (Haf fner et alj|2003 : 
iMadsen fc Revnolddl2005l ). Fi guresQii and [2] show a map 
of EM sin | b\, the emission measure corrected for path 
length through the WIM. 

2.1. Removed H II region sightlines 



i the WHAM Northern Sky Survey for all \b\ > 10° sightlines (panel 
nel b). Maps are centered on I = 120°; dotted lines represent 15° 
dark spots are sightlines contaminated by bright stars which were 

In order to sample only the WIM, we have removed 
from the WHAM survey sightlines that intersect Ha 
emission due to identifiable, discrete structures with 
known sources of ionizing radiation. We identified 17 
large, classical H II regions with |6| > 10° in the WHAM 
survey, then removed sightlines within a circle centered 
on the ionizing star. We chose the radius of the circle 
by eye to encompass the Ha enhancement of the H II re- 
gion. Because many regions are not circular, this crude 
mask removes some sightlines that may properly belong 
to the diffuse WIM. We list the masked H II regions 
in Table [TJ We also removed two regions with signif- 
icantly enhanced Ha emission which are not classical 
H II regions. The Gum Nebula is a large H II region 
centered at I = 258° and b — —5° with a diameter 
of 36°. Most of the Gum Nebula is below declination 
S = -30°, the southern limit of the WHAM survey. 
The Ha emission from the Orion-Eridanus super-bubble 
(1 80° < I < 270°; -50° < b < 0°) is discu s sed in deta il 
bv lRevnolds & Oedenl (fl979h : iHeiles et~aTl (fT999l[2000h : 
and IMadsen et al.l (|2006h . 

The remaining portion of the WHAM survey has 
26, 221 sightlines which sample primarily the WIM — 
plasma with no readily identifiable, discrete source of 
ionizing radiation. A map of these diffuse WIM sight- 
lines is shown in Figure [JJ>. Considerable structure re- 
mains evident, inc lud ing high-latitude fila ments (e. g. 
lHaffner et allll998f ). IMadsen et all (|2006h found that 
the line ratios [N II] /Ha (which is sensitive to the tem- 
perature of the gas) and [S II] /Ha in three discrete 
high-latitude filaments are different from line ratios ob- 
served in classical H II regions and are indistinguishable 
from the ratios in the WIM. Therefore, we consider the 
high-latitude filaments as part of the W IM sample. We 
did n ot remove "WHAM Point Sources" ([Reynolds et alJ 
2005) — 1 pixel enhancements in the WHAM survey, in 
many cases asscoiated with planetary nebulae — because 
the 84 identified point sources contaminate only 91 sight- 
lines and the excess emission measures are much smaller 
than those of the classical H II regions. 

3. DISTRIBUTION OF EMISSION MEASURE 

Histograms of EMsin|6| are shown in Figure [Hh. The 
histogram for the entire high-latitude (|6| > 10°) WHAM 
survey shows a peak near EMsin|6| = 1.4 pc cm -6 . 
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Fig. 3. — Histograms of EMsin |b| for all |6| > 10° sightlines in 
the WHAM survey (thin, solid lines) and for diffuse WIM sightlines 
(thick, solid lines). Gaussian and lognormal fits to the diffuse WIM 
distribution are shown with dashed lines. Neither sample is well fit 
by a normal distribution (panel a), but the WIM is well matched 
by a lognormal distribution (panel b). Residuals comparing each 
fit to the diffuse WIM distribution are shown in the lower panel. 
They are defined as (iVwiM — ^VfitV-^fit, where A^im is the value 
in each histogram bin, are shown below for the WIM distribution. 



TABLE 2 
Lognormal Fits of EM sin 1 61 



range 


^Q(log EMsin 

( pc cm -6 ) 


°"logEMsin \b\ 

(dex) 


N 


x 2 


Full WIM 


1.406±0.004 


0.190±0.001 


26221 


3.9 


10 < |6| < 30 
10 < b < 30 
-30 < 6 < -10 


1.798±0.009 
1.647±0.011 
1.979±0.012 


0.192±0.002 
0.200±0.002 
0.135±0.002 


8747 
5439 
3308 


1.6 
1.3 
4.0 


30 < |6| < 60 
30 < b < 60 
-60 < b < -30 


1.356±0.005 
1.316±0.006 
1.409±0.006 


0.156±0.001 
0.177±0.002 
0.115±0.001 


12296 
8216 
4080 


3.9 
2.9 
1.8 


60 < |6| < 90 
60 < b < 90 

-90 < 6 < -60 


1.068±0.005 
0.991±0.006 
1.198±0.009 


0.139±0.002 
0.132±0.002 
0.127±0.003 


5178 
3246 
1932 


3.0 
3.8 
0.7 



Note. — Gaussian fits to the distribution of log EM sin \b\ for 
subsets of all WIM sightlines - the WHAM survey at \b\ > 10° 
with classical H II regions removed. There are N pointings in 
each region; x 2 is the reduced x 2 goodness-of-fit parameter. See 
Figs. [3] and |1 

When this histogram is plotted with a linear EMsin |6| 
axis, the distribution shows an extended tail to large val- 
ues of EMsin | b\. Restricting the sample to the WIM 
sightlines reduces this high EMsin|6| tail, but does not 
eliminate it. The figure also shows that the EMsin |&| 
histogram is clearly not well fit by a normal (Gaussian) 
distribution; at high EMsin the data can be as much 
as 27 times greater than the best-fit Gaussian. 

However, when the EMsin \ b\ distribution is plotted on 
a logarithmic scale for the horizontal EM sin \b\ axis (Fig- 




FlG. 4. — Histograms of log EM sin 1 6 for WIM sightlines in a 
range of Galactic latitudes with Poisson errors. Lognormal fits 
and residuals (as defined in Fig. [3} are shown; fit parameters are 
shown in Table [2] 

ure[3f)), we find that the distribution is extremely well fit 
by a lognormal distribution of EM sin |6| or, equivalently, 
a normal distribution of logEMsin|6|. This is particu- 
larly true for the pure WIM sample, where the greatest 
deviation between the data and the best fit lognormal 
distribution at high EMsin |6| is a factor of 3.0. In fact, 
the fit is so good that a comparison of a lognormal distri- 
bution to the observed EMsin |6| distribution is a useful 
indicator of the quality of the chosen WIM sample. With 
an incomplete removal of H II regions, the fit to a log- 
normal distribution is poor at high emission measures. 
Parameters of single-component fits of the distribution 
of log EM sin |6| for WIM sightlines are listed in Table El 
The reduced \ 2 goodness-of-fit parameter (e. g. iTaylorl 
[1982), 



x 2 = 



£ 

bins i 



Ni - /(EMsin|6|i) 



/d, 



(2) 



is also listed, where Ni is the number of pointings in 
bin i, (jjvi = \/N~i is the Poisson error, /(EMsin \b\i) is 
the modeled count in each bin, and d is the number of 
degrees of freedom in the fit. We fit Gaussians with a 
zero offset parameter, so each fit has four constraints. 

3.1. Variations with latitude 

We found no evidence for any significant change of 
the EMsin | b\ distribution as a function of Galactic lon- 
gitude, but there is a dependence on Galactic lati- 
tude. This trend is visible in the all-sky map (Fig. []} 
which shows that EMsin |6| appears to be lower at high 
latitudes than at low latitudes. The distributions of 
EMsin | b\ in three latitude bins are plotted in Figure HJ 
lognormal fit parameters are given in Tabled The center 
of the distribution of EMsin |6| is ~ 0.2 dex lower at high 



5 



> 60°) than at low (|6| < 30°) latitude. The distri- 
bution is also ~ 0.05 dex wider at low latitude, although 
this trend is not present at negative latitudes, which are 
less completely sampled by the WHAM Northern Sky 
Survey than positive latitudes. 

In a plane-parallel medium, EMsin|6| would be con- 
stant with latitude. The observed variation with lati- 
tude may indicate that the Sun is in a relatively low- 
density region of the Galaxy because lower-latitude sight- 
lines sample longer pat hlengths within the plane than 
high-latitude sightlines ( Re vnolda 1 1 9 9 7l ) . For example, a 
b = 10° sightline samples a cylinder of radius 5.7 kpc be- 
fore e?sin \b\ is greater than the 1 kpc scale height of the 
WIM, whereas a b = 60° sightline samples only within 
a 0.6 kpc cylinder around the Sun. This low density 
region of ionized gas in the solar environs is probably 
not primarily due to the Local Bubble, which is much 
smaller (~ 100 pc in radius) than the scale height of 
the WIM and thus only accounts for a small fraction 
of the line of sight. Additionally, some lower latitude 
sightlines include spiral arms, which have higher densi- 
ties than the more local gas sampled by high-latitude 
sightlines (jHaffner et all 12003( 1 . 

4. TURBULENCE IN THE WARM IONIZED MEDIUM 

What is the physical mechanism that produces the ob- 
served lognormal distribution of EM sin |6|? If one starts 
with an average density for the WIM and repeatedly 
adds or subtracts a small increment to each volume el- 
ement, the central limit theorem shows that the result- 
ing probability distribution function (PDF) of density 
will be normal (Gaussian). Similarly, a series of small 
compressions and rarefactions, which is expected for a 
compressible fluid, will produce a lognormal distribution 
because multiplicat i on is additive in logarithmic units 
(IVazauez-Semadenil Il994t iPassot fc Vazauez-Semadenil 
Il998t IVazquez-Semadeni fc Passotl Il999h . Thus it is 
quite natural to expect a lognormal density distribution 
in the WIM. Models of turbulence, f or example, produce 
a lognormal distribution of density (iPadoan et all 119971 : 
iNordlund fc PadoaiJl999HOstrike'r et al.lll999Ll2001h . In 
the presence of magnetic pressure, the lognormal distri- 
bution of density is not exact. However, at large Ma, 
the magnetic pressure deco rrelates from density, leaving 
a log normal density PDF (Pass ot fe Vazquez-Semadenil 
1200311 . The effect of magnetic fields on the PDF is also 
small when shocks in the medium are randomly oriented 
with respect to the field (|Beresnvak et al.ll2005( l. 

If the path length of a column is long enough for 
elements within the column to be uncorrelated, the 
central limit theorem implies that the column den- 
sity distribution should approach a Gaussian. How- 
ever, if the turbulence is correlated on a scale longer 
than the path length, the column density is expected 
to have the same distribution fun ction as the density 
(IVazauez-Semadeni fc Garcial |2001( ). as is seen numeri- 
cally ( Ostriker et alj 12001 1. Because log(n 2 ) = 21ogn, 
if the column density (J nds) is lognormally distribued, 
we expect the emission me asure (f n 2 ds) to be lognor- 
mally distributed as well. IVazauez-Semadeni fc Garcial 
(2001) define the ratio of the path length to the decor- 
relation length as the parameter 77, with full conver- 
gence to a Gaussian column density distribution ex- 
pected for 77 1 and partial convergence for 77 > 1. 



The decorrelation length is not know n a priori, although 
IVazauez-Semadeni fc Garcial (|200lH suggest that it may 
be related to the outer scale. 

The outer scale of the turbulence in the WIM in 
the Milky Wa y is not known, with estimates rang- 
ing fr om 4 pc dMinter fc Spanderill996h to 30 - 100 pc 



e. g I Armstrong et al. 1995J ; de*Avillez fc Breitschwerdtj 
120071) . IChepurnovet all (120061) apply the velocity co- 
ordinate spectrum (jLazarian fc Pogosvanl I2006D tech- 
nique to H I data and find an injection scale of 
105 pc. Studies of H I in dwarf irregular galaxies 
such as the Large and Small Magellanic Clouds and 
Holmbcrg II have found outer scales of several kilopar- 
secs, much larger than most c l aims within the Milky 
Way ( Stanimi rovic et all Il999t lElmegreen et al.l 120011 : 
iDib fc Burkertll2005H . iHaverkorn et all (|2006ll2008f) sug- 
gest that the outer scale in the magneto-ionized medium 
in spiral arms is quite small (a few parsecs) with the tur- 
bulence related to H II regions, while the global magneto- 
ionized ISM has a much larger outer scale (~ 100 pc) 
with turbulence driven by supernovae. At a distance of 
1 kpc, the 1° beam of WHAM corresponds to a trans- 
verse length of 17 pc, so our data are not sensitive to 
parsec-scale fluctuations. 

5. MHD MODELS 

We have chosen to compare our observational results 
to models of MHD isothermal turbulence. The WIM is a 
nearly fully ionized plasma and is observed to be magne- 
tized. The isothermal assumption is justified because the 
temperatures of the WIM are observ ed to lie in the range 
6000 - 10000 K (jHaffner et al.lll999l) . so the sound speed 
varies by < 30% within the medium. (|Kowal et al.ll2007l 
hereafter KLB) developed simulations of 3D compress- 
ible, isothermal MHD turbulence over many dynamical 
times for an extended range of sonic and Alfven Mach 
numbers. Although these models were not developed 
with the WIM in mind, they provide a useful context 
for interpreting our observations. 

For a plasma with thermal pressure p, mass density p, 
and magnetic field strength B, each simulation is prin- 
cipally characterized by two numbers: the sonic Mach 
number, M s — v/c s , and the Alfven Mach number, 
Ma — v/ca, where v is the particle speed in each 
cell. The sound speed and Alfven speed are given by 

The 



c s = \J dp I dp and ca = B/^/inp, respectively, 
models are described in detail by KLB and references 
therein. The turbulence was driven using a solcnoidal 
forcing function in Fourier space; this minimizes the im- 
pact of the chosen forcing function on the resulting den- 
sity structures. 

KLB analyzed the statistics of 3D density and 2D col- 
umn density, including PDFs, spectra, structure func- 
tions, and intermittency of density and the logarithm of 
density, in order to establish the relation between the 
statistics of the observables (such as column densities) 
and the underlying 3D statistics of density. They found 
that the amplitude of density fluctuations strongly de- 
pends on the sonic Mach number in both weakly and 
strongly magnetized turbulent plasmas. Moreover, PDFs 
of density fluctuations have a lognormal distribution with 
a width growing with M s - The KLB models have 77 
ranging from 5 (for subsonic models) to 13 (supersonic), 
meaning only partial convergence to a Gaussian column 
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density PDF is expected; the column density distribu- 
tions of the models are very nearly lognormal, like the 
underlying density distributions. 

We match the simulations presented by KLB to the 
physical conditions of the WIM and compare the re- 
sulting distributions of simulated emission measure and 
column density to WHAM and pulsar data. We use 
256 x 256 x 256 element simulations in both super- and 
sub-Alfvenic and super- and sub-sonic regimes. Most of 
the models analyzed here were presented in KLB, but we 
added two mildly supersonic (M s = 1.4—1.7) models for 
this work. Most models were run for at least 10 dynam- 
ical times and all reached approximate steady-states in 
all of the statistical parameters. 

5.1. Setting physical scales 

Since the simulations were performed in dimensionless 
coordinates, it is necessary to use three parameters to 
specify the physical dimensions of the simulation. We de- 
scribe this procedure in detail in Appendix lAl Once the 
physical parameters of each simulation are uniquely spec- 
ified, we can determine how well each simulation matches 
the observed emission measure PDF as well as other ob- 
servations. As described below, we have chosen to fix 
the following three parameters: (1) total vertical disper- 
sion measure through the box, DM sin |6| = 23 pc cm~ 3 , 
(2) length of the box, L box = 200 - 1000 pc, and (3) 
temperature of the gas, T = 8000 K (see §|2J|. 

Dispersion measures of pulsars more than ~ 3 scale 
heights above the midplane sample gas from the entire 
WIM, so we use the value of DM sin |6| for these pulsars 
to set the column density through the simulation cube. 
In practice, the only pulsars high above the midplane 
with reliable distance measurements are those in globular 
clusters. Six gl obular clusters with \z \ > 3 kpc contain 
known pulsars ({Manchester et aLll2005f ) . A recent search 
for pulsars in all globular clusters visible from Arecibo 
failed to detect pulsars in additional clusters at \z\ > 
3 kpc (|Hessels et al.ll2007l ). The mean DM sin |6| for the 
\z\ > 3 kpc globular clusters is 23 pc cm~ 3 (Table [3]). 

Dispersion measures of pulsars with DM-independent 
distance measurements establis h the scale hei g ht of 
the WIM as hm u = 1 kpc (|Revnoldsl I1989L 119971 : 
iGomez et al.l l200lh : Ha i ntensities towards t he Perseus 
arm yield the same result (jHaffner et al.l [T999) . A recent 
analysis of high-latitude (\b\ > 40°) pulsars with distance 
measurements from either parallax or association with a 
globular cluster found a considerably higher s cale height 
for the WIM of 1.8 kpc (|Gaensler et all 12008). Because 
this recent study samples primarily gas within a ~ 2 kpc 
cylinder around the Sun, it may imply a different WIM 
scale height locally than in spiral arms. Comparisons 
of EM and DM have indicated that the WIM occupies 
~ 20 — 40% of the volu me within the 2 kpc thick layer 
about the midplane (e. g. lRevnoldslll99ll but see § 17. ip . 
Combinations of emission measures and dispersion mea- 
sures and dispersion measures and the scale height both 
constrain the number density of the gas, so the problem is 
overconstrained. Therefore, we use simulation box sizes 
ranging from 200 — 1000 pc and test each scale by com- 
paring simulated emission measures to the WHAM data. 
Given that the actual occupied volume of the WIM must 
be no greater than the observed scale height, we do not 
test box sizes > 1.0 kpc even though the scale height 



could be larger. The solenoidal forcing function driving 
the turbulence has a wavelength one-quarter the length 
of the simulation box, which corresponds to 50 — 250 pc. 

5.2. Determining the best-fitting model 

Each of the twelve KLB models used here is uniquely 
identified by two parameters: the initial thermal pres- 
sure, pi n i, and a uniform external magnetic field strength, 
£?ext- In dimensionless units, the values used were -B ext = 
0.1 ("weakly" magnetized) and 1.0 ("strongly" magne- 
tized), and p ini = 0.01, 0.1, 0.2, 1.0, 4.0, and 8.0. Af- 
ter a simulation has reached a steady-state, it can be 
characterized by a final mean sonic Mach number M s 
and mean Alfven Mach number Ma- The sonic Mach 
number decreases with increasing initial pressure; as p; n i 
rises from 0.01 to 1.0, M s drops from 7 to 0.7. The 
Alfven Mach number rises with decreasing B cxt , yielding 
M A = 0.5 - 0.6 for high B oxt cases, and M A = 2.0 - 2.3 
for low £?ext cases. We note that these Mach numbers are 
volume averages. This may overemphasize the impact of 
low density cells which do not contribute strongly to the 
observed properties of the gas but have systematically 
higher velocities than more readily-observed, higher den- 
sity cells. However, the mean Mach number of the dens- 
est 5% of the cells in the simulation cube typically differs 
from the mean over the entire cube by less than 10%. 

Using the physical scaling discussed above, we calcu- 
late simulated emission measures for each model (Fig- 
ure [5|) to determine which model, if any, provides the 
best match to the WHAM data. In each case, the simu- 
lated emission measure is calculated along each of 256 2 
sightlines oriented perpendicular to the external mag- 
netic field. (For the sonic Mach number required to 
match the data, the predicted EM sin \b\ distribution does 
not depend significantly on the orientation with respect 
to the magnetic field, as discussed in § 17.31 ) Parameters 
characterizing the EM sin |6| and density distribution for 
each model are given in Table 01 

Figure [6] shows the mean and standard deviation of 
the lognormal distribution found in the WHAM observa- 
tions compared to Gaussian fits to the lognormal distri- 
bution found in each of the models. A successful model 
must match both the mean, (log EM sin \ b\) , and the stan- 
dard deviation, cri og EMsin |6|, of the WHAM distribution 
of log EM sin The standard deviation is highly sensi- 
tive to the sonic Mach number; higher and wider distri- 
butions of EM correspond to higher sonic Mach numbers. 
Mildly supersonic (1 < M s < 3) simulations best match 
the width of the WHAM distribution. The spread in 
the EM sin | b\ distribution is weakly sensitive to Alfven 
Mach number; models with comparable pressure but Ma 
varying by a factor of four have comparable EM distri- 
butions. Note that the mean magnetic field strength in 
the strongly magnetized, M s ~ 1 — 3 models (column 2 
of Tbl. I} is 2.4 - 3.4 fiG; the WIM i s thought to h ave a 
magnetic field strength of a few /xG (|Ferrierell2001h . 

Once an appropriate value of M s is selected to match 
the width of the log EM sin |6| distribution, the physical 
box size can be varied to match the mean value. For a 
fixed box size, high M s models with a relatively wide dis- 
tribution of EM have higher mean values of EM than low 
M s models. Because EM is the integral of n 2 , the mini- 
mum value of EM for a given column density (J n ds) is 
for a perfectly homogeneous medium, whereas the max- 
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TABLE 3 

Globular cluster pulsar and emission measure data 



Globular 


I 


b 


2 


DMsin|6| 


EM sin |6| 


n c 


L c sin |b| 


Cluster 


n 


n 


(kpc) 


(pc cm -3 ) 


(pc cm -6 ) 


(cm~ 3 ) 


(PC) 


M3 


42.3 


78.7 


10.2 


25.89 ±0.05 


0.95±0.03 


0.037±0.001 


703± 22 


M5 


4.0 


46.8 


5.7 


21.66 ±0.03 


1.31±0.02 


0.060±0.001 


359± 5 


M13 


59.1 


40.9 


4.8 


19.87 ±0.02 


1.12±0.06 


0.056±0.003 


352± 19 


M15 


65.1 


-27.3 


-4.1 


30.812±0.014 


3.54±0.04 


0.115±0.001 


268± 3 


M30 


27.2 


-46.8 


-6.5 


18.283±0.003 


1.94±0.07 


0.106±0.004 


172± 6 


M53 


333.1 


79.8 


16.9 


23.6 ±1.5 


1.14±0.07 


0.048±0.004 


488± 53 


Note. - 


Propert 


les of lines of 


sight towards \z 


| > 3 kpc globular clusters 


with known 



pulsars. EM data from the WHAM survey; sightlines towards M3 and M5 are contam- 
inated by bright stars ( see text). DM and distance data from the ATNF pulsar catalog 
Manchester et al.|[2005T) . 
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Fig. 5. — Histograms of Ha emission measure EMsin|6| from WIM sightlines in the WHAM northern sky survey (thin lines) and 
simulated EM from MHD simulations with a box size of 500 pc (thick lines). Lognormal fits of the simulated EM distribution are also 
shown (dashed lines). Initial external magnetic field strengths and initial pressures computed from each model are specified in code units, 
and sonic and Alfvenic Mach numbers are also listed. 

imum EM for a given column density occurs when the 
gas is heavily clumped. Figure [7] shows the effect of 
changing the box size of the model which best matches 
the observed width of the distribution. A box size of 
Lbcx ~ 400 — 500 pc best matches the mean EM for the 
mildly supersonic models. For the remainder of this pa- 
per, we adopt a box size of 500 pc. 

5.3. Effects of viewing geometry 

In the previous section, we implicitly placed an ob- 
server at different points along the plane at one face of 
the cube to look vertically through the cube with parallel 
sightlines. This visualization, used by KLB, has the ad- 
vantage of simplicity and efficiently utilizes every cell of 
the simulation cube. However, in reality, we only enjoy 



one vantage point, which could potentially be biased by 
nearby features in the WIM. We thus consider the effects 
of a different viewing geometry. 

Here, we place the observer at the center of one face 
of the cube and look through the cube with diverging 
sightlines above a "latitude" of 65°. Sightlines below 63° 
emerge from the side of the cube before sampling its full 
height. This visualization more accurately reflects obser- 
vations of the WIM but is sensitive to the details of the 
vertical distribution of the gas, which is poorly under- 
stood and is not modelled by these simulations. Regions 
with a relatively small volume and an enhanced or de- 
pressed density that are near the observer cover a large 
solid angle and can therefore skew the distribution of in- 



TABLE 4 
Model fit parameters 



Bcxt 
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\ A 
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A4 4 
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i n(los EM) 


""logEM 


,-,2 
X 
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1U > " * 


'-'"log n 


(n c ) 


0"n c 


i J \ 


a L c 




(HG) 








(P C i 


(pc cm - ®) 


{ QCX ) 




( cm - ^) 


(^clcx ) 


( cm - ) 


( cm~^) 


(pc) 


(P c ) 


(1) 


(2) 


(3) 


(A \ 

( 4 ) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 
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Data (Tables |2|&|3): 
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Note. — Parameters of models. (1): External magnetic field. (2): Final mean magnetic field feq. rA8l l. (3): Inital pressure. (4) & 
(5): Final, volume-averaged Mach numbers. (6): Box size. (7) & (8): Most probable EM and standard deviation of log EM. These 
data with a box size of 500 pc are plotted in Figure [6] (9): Reduced x 2 goodness-of-fit parameter of lognormal fit to simulated EM 
distribution. (10) & (11): Most probable density and standard deviation of log n. (12) - (15): Characteristic density and occupation 
length (see §17701. 
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Fig. 6. — Comparsion of observed and modeled EM sin |b| ranges 
as a function of both sonic and Alfven Mach number. We use a 
500 pc box size in each case. Vertical gray bars depict the ±1<t 
range of Gaussian fits to the log (EM sin |6|) distributions of each 
of the twelve models considered. The thinner bars correspond to 
the weakly magnetized (high Ma) cases; the thicker bars cor- 
respond to the strongly magnetized (low Ma) cases. The left 
panel shows both the mean and spread increasing as M a increases. 
The solid and dashed horizontal lines show the observed mean 
((log (EM sin |6| (pc cm - 6 ) -1 )) = 0.146) and standard deviation 
(0.190) of the distribution. The mean of the simulated distribu- 
tions is set by our physical scaling of the simulation (§[5T}, but the 
width is independent of choice of scaling parameters. Error bars 
are omitted for clarity. 
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Fig. 7. — Effects of changing the simulation box size. The 
WHAM WIM sightline histogram is shown with a solid, thick line. 
Three scalings of the Ms =1.7 model are shown (thin lines) with 
box sizes of 200 pc (dashed), 500 pc (solid), and 1000 pc (dotted). 



tensities in this visualization. However, due to the Local 
Bubble, there is little WIM plasma within ~ 100 pc of 
the Sun; placing a gap of 100 pc between the observer 
and the simulation cube lessens the impact of a small, 
nearby clump and yields distributions roughly compara- 
ble to the observed WIM distribution. 

Samples of these visualizations are shown in Figure [8l 
Diverging sightline visualizations yield emission measure 
distributions consistent with those in the parallel sight- 
line visualization. Because the diverging sightline visu- 
alization of the models is limited to \b\ > 65° and re- 
quires detailed knowledge of the vertical distribution of 
the gas to be realistic and because there is no evidence 



that adopting this geometry significantly changes the re- 
sults, we use the parallel sightline visualization for the 
remainder of this paper. 

6. TESTING ADDITIONAL PREDICTIONS OF 
ISOTHERMAL MHD TURBULENCE MODELS 

Although the KLB models of isothermal turbulence 
were not designed for the WIM, we find that they match 
WHAM EMsin|6| data for physically reasonable values, 
namely sonic Mach numbers in the range A4 S = 1.4 — 2.4 
with a box length of Lbox = 500 pc. With the exception 
of the magnetic field strength, which only weakly influ- 
ences the EM sin | b\ distribution, there are no additional 
parameters which are varied in the models used in this 
work. We now present the results of testing model pre- 
dictions of pulsar dispersion measures and Ha velocity 
profiles, which provide additional support for this model. 

6.1. DM distribution 

Pulsar dispersion measures probe the column density 
of free electrons, DM = n e ds, where D is the dis- 
tance to the pulsar (which is known independently of 
DM for relatively few pulsars). Figure [5] shows the model 
DM distributions and the distribution for all pulsars with 
meas ured dispersion measure s in the ATNF Pulsar Cat- 
alog ([Manchester et alj|2005h . obtained on 2006 May 26 
from the pulsar catalog web site. 1 Because pulsars are 
embedded in the medium, DM is a truncated column 
density, while emission measure probes the medium to 
infinity. The scale height of pulsars is poorly known, 
but m ay be compara ble to the 1 kpc scale height of the 
WIM (L orime a fl995h . If there were a sufficiently large 
population of pulsars with |z| > 1 kpc, one would expect 
to see a pile-up of pulsars at the asymptotic value of 
DMsin|6| = 23 pc cm -3 measured with the \z\ > 3 kpc 
globular cluster pulsars. The spread in the asymptotic 
value would be due to the fluctuations in the dispersion 
measure. For the models considered, the predicted total 
column density is reasonably well fit by a lognormal dis- 
tribution with o~ dm, model = 0.03—0.19 over the full pa- 
rameter space. For our preferred models with M. s = 1.4 
to 2.4, we find a\ og dm, model = 0.08 — 0.12. By compar- 
ison, the six globular clusters pulsars are characterized 
by CTlogDMsin |6|,obs = 0.08. 

To test the dispersion measure predictions of the 
model, we randomly placed 10000 pulsars within a simu- 
lated WIM. We assumed that the pulsars have a sech 2 (z) 
distribution (chosen because, unlike an exponential dis- 
tribution, it does not have a cusp at z = 0) and that 
the WIM has a uniform filling fraction up to \z\ = 1 kpc. 
We then placed each pulsar in a random cell in the (two- 
dimensional) modelled column density grid and multi- 
pled the total column density in the cell by z pu i sar //iwiM 
or by one if z pu i S ar > ^wim- Of the 1428 non-globular 
cluster, Galactic pulsars (at all latitudes), 2.4% have 
DM sin |6| > 23 pc cm~ 3 . For each model, we chose the 
pulsar scale height so that the minimum height for which 
the fraction of sightlines with DMsin|6| > 23 pc cm~ 3 
matches this observed value. The required pulsar scale 
height (listed in Figure decreases from ~ 750 pc at 
M s ~ 0.7 to ~ 450 pc at M s » 7 and is ~ 500 pc in our 

1 http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ 
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Fig. 8. — Three visualizations of a strongly magnetized (Ma = 0.6), mildly supersonic (M s = 2.1) model. In panel (c), emission 
measures are shown from the vantage point of each point along one surface of the cube looking vertically in parallel sightlines. In panel (a), 
EM sin |6| is shown in all directions from the vantage point of the center of that plane with a 100 pc-thick plane with density n = below 
the cube. (Sightlines below b = 67°, which do not include the opposite face of the cube, are not shown.) Panel (6) is similar to panel (a), 
but the box has been shifted by half its length in one direction, exploiting the periodic boundary conditions. Corresponding histograms 
are shown below each panel. 

preferred, mildly supersonic models. The required scale 
height is slightly higher in strongly magnetized simula- 
tions than in weakly magnetized ones. 

The resulting distribution of vertical dispersion mea- 
sures is shown in Figure [9] with a thick line. In all mod- 
els, the observed low DM tail and turnover near the ob- 
served peak at DM = 23 pc cm~ 3 are both reasonably 
well matched by the model. Higher M s models under- 
predict the number of DMsin|6| « 20 pc cm~ 3 pulsars 
and over-predict the number of DMsin|6| w 2 pc cm -3 
pulsars. The subsonic and preferred, mildly supersonic 
simulations each produce a DM distribution similar to 
the observed distribution of real pulsars. 

In principle, because we constrain the sonic Mach num- 
ber of the WIM independently of the DM distribution, 
these models allow a measurement of the pulsar scale 
height. However, the crude estimate of the DMsin|6| 
distribution we present here does not consider any selec- 
tion effects in the distribution of known pulsars with dis- 
tance or latitude or multiple disk populations of pulsars. 
Moreover, the resulting pulsar scale height (or heights) 
would likely be substantially different if we employed a 
stratified-density model of the WIM. 

6.2. Velocity profiles 

To determine simulated velocity profiles, we assumed 
a thermal velocity profile in each cell with a center cor- 
responding to the component of the simulated velocity 
along the line of sight. Along a line of sight at coor- 
dinates (a;, y) at a given velocity v, the simulated line 



intensity is 

EM(tf) = J n(x) 2 f(v, x)dz = ^ n(x, y, Zi) 2 f(v, x, y, z t ), 

i 

(3) 

where z is one of the coordinates represented by the po- 
sition vector x = {x,y, z). In a cell at a position x with 
velocity vg, the thermal velocity profile is (eq. 3-20 of 
lSpitzedll978l ) 

-(v~v (S)) 2 /b 2 

f(v,x) = —= , (4) 

0V7T 

where b 2 = iksT/m is the square of the velocity width 
parameter; at 8000 K, b — 11.5 km s - for atomic hydro- 
gen, corresponding to a line full width at half-maximum 
of 19.1 km s~ x . A 1° beam covers approximately 17 pc 
at a distance of 1 kpc, and the size of an individual cell 
in the models is 2.3 pc, so we average over parallel sight- 
lines in a 7 cell square in the simulation cubes to approx- 
imately match a single WHAM beam. We show sample 
simulated line profiles for the mildly supersonic models 
in Figure 1101 Gaussian fits to many line profiles yield 
line widths varying from 20 - 25 km s _1 for M s « 2.4 
models and 18 — 21 km s _1 for M s « 1.4 models. 

How do these values compare to observational con- 
straints? WHAM survey Ha line profiles near the Galac- 
tic pole or towards the Galactic anticenter (I = 180°) 
diagnose the motions of the gas independent of differen- 
tial rotation. The polar regions have coherent regions 
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Fig. 9. — Histograms of pulsar dispersion measures DMsin|f>|, or column densities from the midplane to the pulsar (thin lines). The 
mean DMsin|b| of high-|z| globular clusters, which was used to define the mean column density through the simulations, is shown with a 
dashed line. The results of placing "pulsars" at random heights within the MHD simulations with Lbox = 500 pc and the indicated pulsar 
scale heights (h) are also shown (thick lines; see jj l6.lt . 




Fig. 10. — Simulated line profiles of Ha emission for models 
consistent with the observed EMsin|6| distribution, averaged over 
a 7 cell square within the simulation cube. A single-component 
Gaussian fit is shown for each profile with a dashed line; the full 
width at half maximum of the Gaussian is listed. Dotted lines de- 
pict the thermal line profile at 8000 K (center ed a t v = km s _1 ) 
used to calculate the simulated profiles (see Ej 16.21 1. Intensities are 
in normalized units. 



of rising and falling gas (Fig. 14 of lHaffner etaflT 2003): 
positive velocity (rising) sightlines tend to have lower in- 
tensities than the predominant negative velocity (falling) 
sightlines. In Figure QTJ we show average line profiles for 
four regions. In computing the line profiles, we multi- 
plied the intensity in each pointing by sin |fo| to correct 
for path length through the WIM. For these directions, 
Ha line profiles have a full width at half maximum of 
40 — 50 km s _1 . This is approximately twice as wide as 
the simulated profiles. However, downward moving in- 
termediate velocity gas contributes significantly to the 
observed line profiles, as well as poorly understood re- 
gions of large-scale bulk flows. (The largest scale turbu- 
lent motions would be larger than the box size for our 
simulations, and therefore not predicted by our current 
model.) 

Perhaps a better test of predictions of the model lies 
in choosing lines of sight th at are known to be free 
of such large-scale motions. lHausen et al.l (|2002t ) ob- 
tained Ha line profiles towards two high-latitude sight- 
lines with long WHAM integrations. With 9600 s of 
integration (compared to 30 s integrations used in the 
WHAM survey) towards the O star HD 93521 (I = 
183.1°, b = +62.2°), they identified a local gas compo- 
nent centered at Vlsr = -10 ± 4 km s _1 with a FWHM 
of 22 ±6 km s _1 , where Vlsr is the velocity with respect 
to the local standard of rest. An additional, intermediate 
velocity component centered at Vlsr = —51 ± 3 km s _1 
has a width of 39 ± 7 km s _1 . The sensitive observations 
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Fig. 11. — Ha line profiles from the WHAM survey. Each profile 
is an average of the line profiles (times s in \b \) in the sightlines 
shown in the corresponding map. As in Fig. 1101 a single-component 
Gaussian fit is shown for each profile with a dashed line. Grey lines 
show the spectra used to compute the average line profile. Maps of 
the regions from which spectra were taken are shown at the right. 
Panel (d) is the region known as the Lockman Window. 

allow the identification of separate components which are 
not typically resolved in the lower signal-to-noise survey 
exposures. If the HD 93521 sightline is representative 
and the low- velocity component corresponds to the WIM 
gas, the observed line profiles are therefore comparable 
in width to the simulated profiles. We note that if model 
line profiles had been wider than observations, the model 
would have been definitively ruled out. 

7. IMPLICATIONS FOR THE WARM IONIZED MEDIUM 

7.1. The meaning of "filling fraction" 

Traditionally, density variations in the WIM have been 
characterized by an "occupation length" 



L B = 



DM 
EM 



and "characteristic density" 



EM 
DM 



(5) 



(6) 




of discrete clumps of some fixed density. These param- 
eters assume a bimodal density distribution function in 



n K') 



Fig. 12. — Histograms of density (solid line) and characteristic 
density n c = EM/DM (dotted line) for mildly supersonic models. 

which the density is either n c or zero, with clumps of 
density n c occupying a length L c along the line of sight. 
The filling fraction will depend upon the n(z) distribu- 
tion. For a uniform distribution, the filling fraction is 
/ = L c //iwiMj f° r a plane-parallel, ex ponential distribu- 
tion, / = L c /2hwiM (|Revnoldslll991l ). 

These values for the WIM can be obtained by com- 
bining Ha emission measures and pulsar dispersion mea- 
sures along lines of sight towards pulsars in high-latitude 
globular clusters more than a few scale heights above the 
midplane. For M13, M15, M30, and M53, we used the 
WHAM survey beam which includes the cluster to de- 
termine the EM. Sightlines towards M3 and M5 are con- 
taminated by bright stars (SAO 82944 and SAO 120946, 
respectively, which each have magnitudes V < 7); Fraun- 
hofer lines from the stars cause an underestimate of the 
emission measure along those sightlines in the WHAM 
survey. Therefore, we used an average of sightlines with 
beam centers < 1.5° from the globular clusters, exclud- 
ing the central beam, to estimate the EM. For the six 
globular clusters, the mean values are (n c ) ~ 0.07 cm -3 
with a filling fraction of rj 30% of a uniform density WIM 
disk (14% of an exponential disk) w ith a scale hei ght of 
1 kpc, consistent with the results of lRevnoIdi (|1991h . 

Consideration of a lognormal density distribution mod- 
ifies this physical interpretation. The filling fraction 
within the simulation cube is one, but the density within 
the cube has a lognormal distribution of values rather 
than a single value. The total filling fraction of the WIM 
within the Galaxy is now the volume which the simula- 
tion cube must occupy in order to reproduce the observed 
emission measure distribution divided by the total vol- 
ume of the Galaxy within the layer in which the WIM re- 
sides. In place of the occupation length, L c , and charac- 
teristic density, n c , we have the simulation length, Lbox, 
and a distribution of densities characterized by the most 
probable density, n p k = 10^ og "\ and the standard devi- 
ation, CTio gn (columns 10 and 11 of Tbl.|U respectively). 
Using the synthetic EM and DM values derived from the 
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models, we also calculate synthetic n c and L c values. 
Over the mildly sub/super-sonic range, M s = 0.7 to 2.5, 
n p k/n c ranges from 0.95 to 0.4, and L c /L\> OK ranges from 
0.95 to 0.62. (For the highly supersonic case, these ra- 
tios are ~ 0.1 and 0.4, respectively.) These changes are 
easily understood. Figure [12] shows histograms of the 
characteristic density for 256 2 independent lines of sight 
through the mildly supersonic models. The characteris- 
tic density is considerably higher than the most probable 
density; n c is weighted to higher densities because of the 
n 2 dependence of EM. 

The effect of the different density distribution is to 
increase the volume that one would say is occupied by 
the WIM. Using the best fit box size, L DO x = 500 pc, 
and the observed WIM scale height, /iwiM = 1000 pc, 
50% of the volume within the 1 kpc thick layer above 
the midplane is occupied by the WIM with a most likely 
density of 0.03 cm" 3 . If the t rue scale height is /iwiM = 
1.8 kpc (|Gaensler et al.ll2008f) . the filling fraction would 
be smaller but our study's requirement for the box size 
would not change. In an exponentially stratified model, 
we expect that the mean density at the midplane would 
be similar to that in the uniform density model used here, 
but the filling fraction would be approximately half this 
value. 

In reality, the WIM consists of an unknown number 
of pockets embedded in an ISM with cold and warm 
neutral and hot ionized regions, rather than the single, 
isothermal box employed here. Although the WIM itself 
is observationally approximately isothermal (§ [2]) , inter- 
faces and mixing with other phases may affect the gas 
in way s not modelled here. Ide Avillez &: Breitschwerdtl 
(2004) model the Galactic ISM with supernova explo- 
sions energizing a four-phase medium with heating and 
radiative cooling and find a combined volume filling fac- 
tor of - 0.6 for the "cool" (10 3 < T < 10 4 ) and "warm" 
(10 4 < T < 10 5 - 5 ) gas. They find that pressure equilib- 
rium between the phases does not hold; instead, turbu- 
lent mixing leads to a dynamical equilibrium. There is 
also a large range of temperatures that result, includ- 
ing gas that should be thermally unstable. However, 
this and similar works do not generally take into account 
the ionization state of the gas by tracking the collisional 
and photo-ionization, recombination and radiative trans- 
fer through the gas. Presumably the gas identified as 
"cool" and "warm" in these simulations are a mixture of 
warm neutral and warm ionized components. Thus it is 
difficult to know what the density distribution of the ob- 
servationally selected warm ionized gas would be in such 
a model. 

7.2. Power requirement for a turbulent warm ionized 

medium 

All of the KLB simulations we have considered were 
driven with the same forcing function. After the simula- 
tions reach a statistical steady-state, the energy injection 
equals the energy "radiated" away in the simulation. For 
isothermal simulations, the energy loss rate is assumed 
to be narrowly localized in the region of shocks and com- 
pressions rather than determined by the balance of heat- 
ing and cooling. All of the simulations were driven with 
a power (in dimensionless units) of unity. As shown in 
the Appendix, this can be converted to a physical value. 

For our best-fitting simulations (A4 S = 1.4 — 2.4, 
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Fig. 13. — Effect of orientation of line of sight with respect to the 
magnetic field. Thin lines show the observed, WHAM distribution. 
Thick lines show the modeled EM distribution perpendicular to 
(solid lines) and parallel to the magnetic field (dotted and dashed 
lines) . 

with Lbox = 500 pc), the gas reaches a turbulent 
steady-state with a kinetic and magnetic energy den- 
sity of Ukin = (1.5 — 3) x 10~ 13 ergs cm -3 and U mag — 
(0.8 — 1.5) x 10~ 13 ergs cm" 3 . The power density 
necessary to provide this turbulence is P = (1.8 — 
5.2) x 10~ 25 ergs s _1 cm -3 . Multiplying this value by 
Lbox yields the power per area per side of the plane, 
(3 — 8) x 10 -4 ergs s _1 cm" 2 . This is approximately 
the same power provided by 10 51 ergs of supernova en- 
ergy per 50 years spread over a disk of radius of 15 kpc 
(3 x 10~ 4 ergs s _1 cm" 2 ), but perhaps a bit on the high 
side. However, given the uncertainties in the spatial 
distribution of supernova energy input and the amount 
of energy converted into kinetic energy, these estimates 
should be viewed cautiously. 

Given the limitation of the isothermal assumption in 
our models, do our values provide upper or lower limits 
on the true power requirement of the WIM? The answer 
to this is not clear. One parameter in our models that 
would affect the power requirement is the driving scale, 
which was a quarter of the simulation box size. Using a 
smaller driving scale would probably increase the power 
requirement to maintain the same level of turbulence. 
Replacing the isothermal assumption (which assumes gas 
radiates immediately after shocks) with both supernova 
and photoionization heating and radiative cooling could 
raise or lower the power requirement for the WIM de- 
pending upon the temperature and density structure of 
the resulting medium. For the purpose of this work, it is 
comforting that the power requirement is not excessive. 

7.3. Magnetic field structure of the warm ionized 
medium 

Figure [13] shows EM PDFs for each of the three possi- 
ble orientations of the box (two perpendicular to the field 
and one along the field) for a weakly and strongly mag- 
netized, mildly supersonic model. For each model, the 
three PDFs are similar, with counts slightly higher in 
the parallel view than in the perpendicular at low emis- 
sion measures in the weakly magnetized model. 

Because the distribution of density, column density, 
and emission measure in the models is relatively insensi- 
tive to the magnetic field strengths we considered, these 
data alone cannot effectively constrain the properties of 
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the magnetic field. However, our study's density struc- 
ture of the WIM effectively fixes the density, one of the 
free parameters in interpreting Faraday rotation studies 
of the ionized medium, which depends on the convolu- 
tion of both the density and the magnet ic field structure 
of th e ISM at high Ga lactic latitude (iBoulares fc Coxl 
1990t iZweibel fc Heilesl fl997l iFerrierd [20011: iHanl 120061 : 
de Bruvn et al.ll2006[ ). Comparison of MHD models such 
as those we considered to the combination of density and 
magnetic field probes will be a powerful way to study the 
magnetic field structure of the WIM in the local solar 
neighborhood. 

We have checked to make sure that the magnetic fields 
used in our models included reasonable magnetic field 
stengths after normalizing the models to the data as de- 
scribed above. For the best fitting cases, the mean mag- 
netic field strengths of the models considered are between 
0.8 — 3.4/iG, similar to observational results (|Ferrierd 
2001). A more extensive grid of magnetic field strengths 
would be desirable for comparison with Faraday rota- 
tion measure studies of high latitude extragalactic radio 
sources. 

8. SUMMARY AND FUTURE DIRECTIONS 

Using data from the WHAM Northern Sky Survey with 
sightlines that intersect classical H II regions removed, we 
find that the Ha emission from the WIM has a lognormal 
distribution of EMsin|6|. Fits describing the observed 
distribution are listed in Table [2l If the WIM were a 
uniform, plane-parallel layer, EM sin |6| would be uniform 
at all latitudes, but the emission at \b\ > 60° is roughly 
0.14 dex lower than the emission in the range 10 < \b\ < 
30°. 

Because the WIM is an isothermal, magnetized plasma 
and a lognormal distribution arises from a series of com- 
pressions and rarefactions in a fluid, MHD models of 
isothermal turbulence appear to be applicable to the 
WIM. KLB recently completed such a simulation, which 
we apply to the physical scales and conditions of the 



WIM. We find that the model effectively reproduces the 
observed lognormal distribution of Ha emission mea- 
sures with mildly supersonic (Ai s ~ 1.4 — 2.4) turbu- 
lence, provided that the turbulent plasma fills roughly 
half of the 2 kpc-thick WIM layer. The model also suc- 
cessfully matches constraints from pulsar dispersion mea- 
sures and Ha line profiles. The model results in the WIM 
plasma filling roughly twice the fraction of the volume of 
the layer found in previous, relatively simplistic models 
(|Revnoldsl 1 1 99 lh . with the most probable local density 
(~ 0.03 cm" 3 ) roughly half that found in previous mod- 
els. We explore some of the ramifications of this result, 
including a discussion on the meaning of filling factor, 
the power requirement of the turbulent WIM, the driv- 
ing scale of turbulence, and implications for studies of 
the magnetized WIM. 

The KLB model used here was not designed specifically 
to describe the WIM; it is a generic model of isothermal 
turbulence. The model accurately reproduces the obser- 
vations we have tested, but the solution is not unique and 
other models should certainly be tested against our ob- 
servational result. Probably the most important change 
to be considered in the future would be the develop- 
ment of vertically stratified models. Replacement of the 
isothermal assumption (appropriate for the WIM, which 
is observed to be nearly isothermal, but not the general 
ISM) with supernova and photoionization heating and 
radiative cooling may also be instructive. 



We thank the referee, Enrique Vazquez- Semadeni, for 
a review which led to substantial improvements in the pa- 
per. WHAM is supported by the National Science Foun- 
dation through grants AST 0204973 and AST 0607512. 
G. K. and A. L. are supported by NSF grant AST 
0307869 and the Center for Magnetic Self-Organization 
in Astrophysical and Laboratory Plasmas. 

Facility: WHAM 



APPENDIX 
PHYSICAL SCALE OF MODELS 

The simulations used in this work are conducted in dimensionless units. In order to compare observations directly 
with the simulations, we must impose reasonable physical constraints on the simulations. We adopt the notation that 
a tilde above a quantity denotes the quantity in dimensionless (simulation) units and a subscript zero denotes the scale 
factor; the symbol with neither a zero nor a tilde is the quantity in physical units. A modelled quantity in physical 
units is a = a$a. Observed constraints are denoted with the subscript fobs', and initial conditions that vary with the 
models are denoted with the subscript 'ini'. 

The simulations used in this work output a density p and x, y, and z components of the magnetic field B XtV ^ z and 
velocity v X: y jZ for each cell with grid coordinates (x, y, z). We specify, from observations, three physical scales: a length 
xo , a velocity vq , and a mass density po ■ These three scales are set by specifying the physical size of the simulation 
volume, the isothermal sound speed of the gas, and the dispersion measure of ionized gas through the simulation 
volume, respectively. 

Length scale 

Since each grid cell is a cube, the volume of a cell is 1 unit 3 ; the total volume of the simulation in dimensionless 
units is therefore N x x N y x N z , where for the simulations considered here N x — N y = N z . By specifying a physical 
length of the box, Lbox, we write the length scaling factor 



Note that we treat the box size as a free parameter because the problem is overconstrained and the filling fraction 
of the WIM is not well known. The chosen box size must be smaller than the observed height of the WIM to be 
physically reasonable. 
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Velocity scale 

The velocity scale factor comes from relating the sound speed in code units to that in physical units, 



Cs,obs / Pini / a r>\ 

Cs y Pini 




where p is the thermal pressure. 
The physical sound speed of an ideal gas is 



where T Q b s is the temperature and p is the mean particle mass. The numerical value comes from adopting m = 
2.36 x 1CT 24 g as the mean mass per hydrogen atom, appropriate for a solar metallicity gas, and p = m/2y e , where 
y e = n e /riH = 1.1 for a plasma in which the hydrogen is fully ionized and helium is singly ionized. This yields 
c s = 10.15 km s" 1 (T/8000 K) 1 / 2 . 
For the simulation, the dimensionless sound speed of each model is 



Initially, all cells are set to have a density p ln i = 1; because there is no mass flow in or out of the simulation cube, 
the mean density remains unity at all times, although we leave it in the equations for dimensional clarity. The initial 
pressure p- ln i is an input parameter for each model. The pressure is calculated in each cell by an isothermal equation of 
state, so pj p is constant over time across the entire grid for a given simulation, and the sound speed is c s — (pini/pini) 1 ^ 2 
everywhere. 

Density scale 

The density scale factor comes from relating the surface mass density in the simulation to the surface mass density 
inferred from the observed dispersion measure and is given by 

m(DMsin|6|)obs 

Po = t : • (A3) 

2/e^obsPini 

Since each cell has a unit length in dimensionless coordinates, the dimensionless surface mass density through the 
simulation is 

S = ^2 Pi = piniN x . 

i=0 

The physical surface mass density can be derived by comparison to the mean vertical dispersion measure through the 
half-thickness of the warm ionized medium, where we use the dispersion measure (DM sin |6|) bs of pulsars in globular 
clusters more than 3 kpc above the midplane. This corresponds to a surface mass density of 

m(DMsin|6|) obs 
-obs = / pds= n H mds= — 



= pals = I 
Jo Jo 



Ve 



We can relate the dimensionless and physical surface mass density to get the surface density scale factor, S = S Q b s /S. 
Using the fact that So = Po x o, we get the mass density scaling factor as above. 

Additional scale factors 

The three scale factors above allow us to determine the values of the scale factor for other physical values including 
time, 



t = X ^ = -^J^ (A4) 

V iVxC s ,obs V Pi- 



'-in-} 



pressure and energy density (which have the same units), 

D] 

KhsVeP 



m(DMsin|6|) ob s< obs 
Po = Povo = 7 ? (A5) 



power density, 



Pa m(DMsin|6|) obs c3 hs N x 
Po = — = r — — , (A6) 

to J/e^obsPiniPini 
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and magnetic field strength 

o / 2 /m(DMsin|6|) obB 

S = VPoWo = c s ,obsW r ■ (A7) 

Although we do not explore the effect of magnetic field strengths in this paper (because the density PDF is fairly 
insensitive to this parameter), we note that the mean magnetic field strength that develops in these simulations can 
be obtained by comparing the sonic and Alfven Mach number recorded for each simulation and is given by 

(B) = 0.519 ( g.) (—^) V2 ( — MG, (A8) 
V ' \MaJ VO.Ol cm- 3 / VlOkms" 1 / V ' 

where n p k is the most probable electron density. This field is a mixture of the originally imposed uniform component 
and the random component that developed as a result of the turbulent driving. 
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